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Quantum  equations  of  motion  describing  the  energy  trans¬ 
fer  dynamics  via  intramolecular  anharmonic  coupling  are  pre¬ 
sented  and  solved  numerically.  Dynamical  features  of  the 
average  excitations  of  a  homogeneous  system  (SFg)  and  a  hetero¬ 
geneous  system  (H/W)  are  quantitatively  discussed  in  terms  of 
the  intramolecular  vibrational  relaxation  (IVR)  rates.  Possi¬ 
ble  mechanisms  of  surface-enhanced  desorption  via  low-power 
laser  radiation  are  proposed. 


1.  Introduction 


The  selective  nature  of  laser-induced  processes,  which 

crucially  depend  on  the  intramolecular  vibrational  relaxation  (IVR) 

rate  and  on  the  pumping  rate  of  the  laser  radiation,  have  been 
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widely  studied  for  homogeneous  (gas-phase)  systems  and  more 
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recently  for  heterogeneous  systems . 

For  homogeneous  systems,  several  types  of  selectivity  have 
been  discussed  in  terms  of  the  relative  magnitudes  of  the  pumping 
rate,  the  IVR  rate,  the  intermolecular  vibration-vibration  and 
vibration-translation  rates:  (1)  mode  selective  (for  fast  pumping), 
(2)  molecule-selective  (for  fast  IVR  rate) ,  (3)  nonselective 
vibrational  excitation  (for  strong  intermolecular  vibration-vibra¬ 
tion  coupling)  and  (4)  thermal  excitation  of  the  system.1 

For  heterogeneous  systems,  e.g.,  molecules  (atoms)  adsorbed 
on  a  solid  surface,  the  nature  of  the  laser  selective  effects  are 
characterized  not  only  by  the  intra-  and  intermolecular  relaxation 
processes  of  the  adsorbed  species  and  the  pumping  rate  of  the 
laser  field,  but  also  by  the  surface-induced  damping  factor  (via 

multiphonon  processes)  and  the  laser-induced  surface  dynamical 

8-n 

phenomena,  e.g.,  migration,  desorption  and  predissociation. 

The  laser-stimulated  surface  processes  then  may  be:  (1)  mode 
selective,  (2)  adspecies  selective,  (3)  migration  (lattice-site) 

8 

selective,  (4)  local  heating  and  (5)  nonselective  thermal  heating. 

In  a  homogeneous  system,  the  mode-  and/or  molecule-selec¬ 
tive  excitations  are  possible  by  a  high-power  short  pulse  when 

1  12 

the  IVR  rate  is  small.  '  For  a  heterogeneous  system,  the  laser- 
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stimulated  selective  surface  bond  breaking  (desorption)  is  also 

possible  when  the  phonon- induced  damping  rate  is  small  for  high- 

5  13 

order  multiphonon  processes.  '  The  IV R  rate  therefore  plays 
the  essential  role  in  the  nature  of  the  laser-selective  excita¬ 
tions,  both  for  homogeneous  and  heterogeneous  systems.  In  this 
Letter,  the  dynamical  features  of  the  intramolecular  modes 
governed  by  the  intramolecular  anharmonic  couplings  are  studied 
for  two  systems:  (A)  SFg  (gas  phase)  and  (B)  H/W  (hydrogen  atom 
adsorbed  on  a  tungsten  surface) . 

2.  Quantum  Equations  of  Motion  for  a  Nonlinear  Open  System 

The  total  Hamiltonian  describing  an  anharmonic  quantum 
oscillator  subject  to  laser  (IR)  radiation  of  an  open  system 
may  be  written  as 

pf  f 

H(t,Q1,Q2,««*Q1J)  =  Hq  +  XHanh  +  H'(t),  (1) 

pf  £ 

where  HQ  is  the  effective  unperturbed  Hamiltonian  of  the  system 
(with  N  normal  modes)  ,  *Hanil  is  the  intramolecular  anharmonic 
coupling  and  H'(t)  is  the  laser  field  interaction  Hamiltonian. 

In  a  second-quantization  representation,  the  individual  terms  of 
the  total  Hamiltonian  may  be  expressed  by 
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In  eq.  (2. a),  the  imaginary  part  of  the  complex  effective  Hamil¬ 
tonian  T j  takes  account  of  the  level  broadening  of  the  open  sys¬ 
tem  due  to  factors  such  as  collisional  (pressure)  broadening  for 
a  homogeneous  system  or  surface-phonon-induced  damping  for  a 
heterogeneous  system.  In  eq.  (2.b),  two  typical  terms  of  the 
anharmonic  coupling  are  shown:  the  first  term  representing 
T2  phase  relaxation  processes,  where  no  vibrational  energy  among 
the  intramolecular  modes  is  exchanged,  and  the  second  term 
representing  the  energy  relaxation  processes  where  the  photon 
energy  is  deposited  in  the  active  mode(s)  and  populated  in  all 
the  other  inactive  mode  via  anharmonic  coupling.  In  eq.  (2.c) , 

Vj  (t)  «Vj0  -cos  uit  is  proportional  to  the  derivative  of  the  dipole 
moment  and  the  effective  local  electric  field  of  the  j-th  mode 
with  laser  frequency  w.  We  note  that  the  anharmonic  coupling 
factors  and  are  related  to  the  p-th  derivatives  of  the 

anharmonic  potential  energy  (with  respect  to  the  normal  coordi¬ 
nates)  with  p  =  £+m+n  being  the  order  of  the  multiquantum  coupling 
In  order  to  study  the  energy  transfer  among  the  intra¬ 
molecular  modes,  we  shall  first  set  up  the  equations  of  motion,  of 
the  energy-related  relevant  quantities .  The  solution  of  the  time- 
dependent  Schrodinger  equation 


with  the  total  Hamiltonian  H  given  by  eq. 


(1) 


may  be  written  as 
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where  A(t)  and  B(t)  are  time-dependent  coefficients. 


Substituting  eq.  (4)  into  eq.  (3)  with  the  expression  of  the  total 
Hamiltonian  [eqs.  (1)  and  (2)],  we  find  the  equations  of  motion  of 

N)  to  be 

A,  ****■*, 


the  amplitude  functions  Bj(j=l,2, 


^  **  -*f  (i  Azt)  l 


+  (W)*rF4‘*)-uriA)Bi  , 


(S) 


where  the  external  detuning  (A^)  and  the  internal  detuning  (Aj)  are 

defined  by  Ajj,=  2tt  (v  .  .-v.)-oj  and  AT=2ir  (Jlv  .+mv  .-nv  )#  In  deriving  eq.  (5)  , 

we  have  used  the  usual  rotating-wave  approximation  and  operator 
14 

algebra  such  as 


(6*a)]  p U)  [4*f  ]  ~  F (tX+&-  ^ ) 


where  F (a  )  is  any  operator  function  and  B ^  are  c-numbers . 

The  important  features  of  the  quantum  equations  of  motion 
[eq.  (5)]  are:  (1)  the  j-th  normal  mode  of  the  system  is  strongly 
coupled  to  the  laser  field  when  VjQ+0  and  the  external  detuning  is 
not  far-off  the  resonance  condition  aE=0  (due  to  the  anharmonicity , 

Cj 

we  shall  find  later  that  the  optimal  detuning  is  red-shifted,  i.e., 
Ag>0);  (2)  the  intramolecular  coupling  is  governed  by  the  coupling 
strengths  *G£mn  an<*  t*le  internal  detuning  A^.:  for  the  intra¬ 

molecular  coupling  to  be  significant,  we  shall  expect  a  near 
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resonance  A^O,  which  was  also  shown  by  a  classical  treatment;  ' 

(3)  by  eq.  (5)  and  its  complex  conjugate,  we  find  that  the  equation 
of  motion  for  BjBj  contains  no  T2  dephasing  term,  i.e.,  no  vibra¬ 
tional  energy  is  exchanged  among  the  intramolecular  modes  due  to 
the  T2  dephasing  processes;  (4)  the  amplitude  function  is  related 
to  the  average  excitation  in  the  Poisson  distribution1® 

j 

where  Pn  is  the  energy  population  of  a  harmonic  oscillator  (for 
A=0)  and  n^  and  are  the  vibrational  quantum  number  and  the 
average  excitation  of  the  j-th  mode,  respectively;  and 

^  2 

Therefore  the  average  excitation  n | B ^ |  is  a  relevant  quantity 
in  describing  the  energy  distribution  among  the  intramolecular 
modes . 

3 .  Dynamical  Excitation  and  Optimal  Detuning 
(A)  SFg  gas-phase  system 

To  show  the  role  of  the  external  and  internal  detunings  in 
IVR,  we  now  first  consider  the  gas-phase  SFg  system  in  the  colli¬ 
sionless  low-pressure  regime  (i.e.  1^=0)  and  focus  on  the  case  of 

2 

the  fourth-order  anharmonic  coupling  EQ3Q2  Qg»  w^ere  e  1S  t*ie  an_ 
harmonicity  and  Q^t(^2  an<*  are  nonna^  coordinates  of  the 
(active),  and  vg  modes.  The  corresponding  internal  detuning 
then  is  AI=2ir  (v2-2v2“Vg)  and  the  external  detuning  is  AE=2irv.j-w» 
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where  ui  is  the  frequency  of  the  tunable  C02  TEA  laser  which  is 
near  resonant  to  the  mode.  The  quantum  coupled  equations  of 
motion  [eq.  (7)]  for  this  case  become 


jts,  =  M 

it&r  Xc+&*t£+2>!ctc*Q a, g&j  r9ej 

.<t  4  =  x '.Q*gg+xc,  S9.CJ 


where  Bj(j=l,2,3)  are  the  related  amplitude  functions  of  the 

normal  modes  v3,  v2  and  vg,  respectively;  X'=12X,  A^=2n (v^  +  vg-2v2) 

and  A '=27rvl-aj,  with  v!(j=2,3,6)  are  the  anharmonic-corrected 
E  j  3 

frequencies,  i.e.,  v ,!=v . -X 'C ^/H;  are  the  quantization  con- 

1/2 

stants  C  .  =  (1f/2m  .a)  ■ )  .  The  above  system,  for  the  case  of  exact 

•J  J  J 

resonance  A'  =0,  was  discussed  and  the  related  IVR  was  estimated 

L* 

12 

by  Gan  et  al.  Here  we  shall  show  the  dynamical  features  of  IVR 

and  find  the  optimal  external  and  internal  detunings  by  solving 

the  coupled  equations  of  motion  numerically.  We  investigate  the 

time  evolution  of  the  average  excitations,  n1=  |  B1 1  ,  xi2=  |  B2  |  and 
_  2 

H3=|b3|  ,  governed  by  eqs .  (9)  with  the  initial  Boltzmann  popula¬ 
tions  at  room  temperature  (T=300K) ,  B^(0)=0.10,  B2(0)=0.214  and 
B3(0)=0  417,  and  the  initial  phases  tan  ^ [I j (0)/Rj (0) ]=0,  where 
1^  and  Rj  are  the  imaginary  and  real  parts  of  Bj(j=l,2,3). 

Fig.  1  shows  the  peak  value  of  the  average  excitation 
(ix*)  of  the  active  (v3)  mode  as  a  function  of  the  external  detun- 


8  2 

ing  with  laser  intensity  1=10  W/cm  .  It  is  seen  that  the  optimal 

external  detuning  A*=  17  cm-1  and  the  value  of  n*  decreases  rapidly 

»  * 
tor  A£  larger  than  A£.  This 


1 
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"red-shift"  behavior  of  the  optimal  detuning  is  a  general  feature 
of  any  laser-excited  anharmonic  oscillator  where  the  energy  spac- 
ings  are  closer  for  higher  levels.  [For  more  detailed  analysis 

we  refer  the  reader  to  Ref.  17.]  Figs.  2(A)  and  (B) 

show  the  exact  resonance  and  the  far-off  resonance,  respectively. 

It  is  seen  that  the  peak  value  of  the  v2  mode  (ji£)  decreases  from 
n*~0.4  (for  A^O)  to  n**0.1  (for  A^=63  cm  1)  .  By  our  numerical  results, 
we  estimate  that  the  rising  times  of  the  average  excitations  tj 
[defined  by  the  first  peak  values  of  the  average  excitation  nj(t*)] 
range  from  t*=27  picoseconds  (ps)  to  t£=70  ps  depending  on  the 
internal  detuning  A^.  However,  the  rising  times  of  the  v2  mode, 
i.e.,  the  IVR  times  are  in  the  order  of  t*  =  10  ps  (or  3.3  cm-1). 

A  typical  excitation  diagram  is  shown  in  Fig.  2(C)  for  the  case 
of  (A^/27T , Aj/2tt)  =  (2.71,1.00). 

Finally,  we  consider  the  optimal  internal  detuning  where  the  intra¬ 
molecular  anharmonic  coupling  is  the  most  significant  (or  the  maxi¬ 
mum  excitation  of  the  v2  mode  n|) •  We  see,  from  Fig.  3,  that  the 
optimal  internal  detuning  is  A j/ 2^*1. 6  cm  1  when  the  external 
detuning  is  fixed  at  its  optimal  value  A/e/2tt,“2.7  cm  . 


(B)  H/W  heterogeneous  system 

~1  1 8 

The  IR  spectrum  of  H/W  is  given  by  (in  units  of  cm  )  : 


^=1048. 6,  v2=645.3  and  v3=1290.6.  The  fourth-order  anharmonic 
coupling  is  given  by  e'Q^Q2Q3,  i.e.,  Ai=2tt  (2v^-v2~v3)  ,  where  the 

mode  is  vibrationally  excited  by  a  C02  laser.  Again,  we  investi¬ 
gate  the  average  excitations  of  these  three  modes  for  the  case  of 


Tj-O.  Since  the  surface-phonon-induced  level  width  (Tj)  simply 
causes  the  decay  of  the  energy  excitations,  it  does  not  affect 


the  rates  of  IVR  significantly.^ 
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Two  typical  dynamical  features  of  the  average  excitations 

n^/  n2  and  JI3  (for  v^,  v2  and  v3  modes,  respectively)  with  laser 

8  2  "1 
intensity  1=10  W/cm  and  the  anharmonocity  e'=5  cm  are  shown 

in  Figs.  4(A)  and  (B) .  It  is  seen  that  these  dynamical  patterns 

of  the  excitations  governed  by  e'Q^^Qj  are  signif icantly  different 

from  that  of  the  SFg  system  governed  by  eQ3Q2Qg.  We  note  that 

•y 

the  former  system  involves  the  operator  function  aia2a3  corres¬ 
ponding  to  a  two-photon  processes  followed  by  an  anharmonic  two- 
quanta  coupling,  while  the  latter  system  involves  a3(a+)  a+  cor¬ 
responding  to  a  single-photon  processes  followed  by  an  anharmonic 
three-quanta  coupling.  We  also  note  that  the  IVR  rate  of  the 
H/W  system  is  faster  than  that  of  the  SFg  system;  however,  the 

average  excitation  of  the  active  mode  in  H/W  system  is  lower  since 

2 

the  two-photon  transition  (governed  by  y^EQ  )  is  usually  less 
probable  than  the  single-photon  transition,  i.e.,.  the  dipole 
transition  governed  by  y^EQ,  where  y^  and  yJJ  are,  respectively, 
the  first  and  second  derivatives  of  the  active  dipole  moment 
with  respect  to  the  normal  coordinate  Q  and  E  is  the  local  elec¬ 
tric  field  acting  on  the  species. 

In  conclusion,  we  note  that  for  a  heterogeneous  system 

the  local  electric  field  acting  on  the  adspecies  may  be  greatly 

19 

enhanced  by  the  surface  effects  such  as  roughness.  Furthermore, 

the  surface  bond  breaking  can  also  be  possible  via  van  der  Waals 

20  21 

interaction  or  the  surface-migration-induced  ra£e  processes. 

One  of  us  (J.L.)  would  like  to  thank  Prof.  Z.  Gan  of 
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Figure  Captions 

Fig.  1.  The  peak  value  of  the  average  excitation  of  the  active 

mode  (n*)  vs.  the  external  detuning  (A^) .  [SFg  system]. 

Fig.  2.  The  time -dependent  average  excitations  of  the  active 
mode  (n1)  ,  the  v2  mode  (n2)  an d  the  vg  mode  (H-j)  for 
the  detunings  (A^/2tt  ,  Aj/210  =  (A)  (0,0)  exact  resonance, 

(B)  (0,20)  far-off  internal  resonance  and.  (C)  (2.71,1.00) 

A7 /2ir  =  2.71  is  the  optimal  external  detuning. 

Fig.  3.  The  peak  value  of  the  average  excitation  of  the  v2 

mode  (n£)  vs.  the  internal  detuning  (Aj)  with  the  opti¬ 
mal  external  detuning  AE=A*. 

Fig.  4.  The  dynamical  features  of  the  average  excitations  of 

the  active  mode  (n^)  »  the  v2  mo<^e  ^2^  an<^  V3  mo<^e 
(fJ3)  for  a  heterogeneous  system  H/W  with  detunings 
(A7 /2tt,A/t/2tt)  =  (A)  (0,0)  exact  resonance  and  (B) 

(10,80) . 
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